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Measuring the distribution of current fluctuations through a 
6c Josephson junction with very short current pulses 
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Abstract. - We propose to probe the distribution of current fluctuations by means of the escape 
probability histogram of a Josephson junction (JJ), obtained using very short bias current pulses 
in the adiabatic regime, where the low-frequency component of the current fluctuations plays a 
crucial role. We analyze the effect of the third cumulant on the histogram in the small skewness 
limit, and address two concrete examples assuming realistic parameters for the JJ. In the first one 
we study the effects due to fluctuations produced by a tunnel junction, finding that the signature 
of higher cumulants can be detected by taking the derivative of the escape probability with respect 
to current. In such a realistic situation, though, the determination of the whole distribution of 
current fluctuations requires an amplification of the cumulants. As a second example we consider 
magnetic flux fluctuations acting on a SQUID produced by a random telegraph source of noise. 



Introduction. — The electronic transport properties 
of a mesoscopic conductor are completely characterized by 
its Full Counting Statistics (FCS), introduced in Ref. [1] 
and defined as the probability distribution for the transfer 
of charges over a certain time interval. At zero tempera- 
ture and in the absence of interaction, for example, quan- 
tum transport is characterized by a binomial distribution. 
This contrasts with the Poissonian statistics, characteris- 
tic of classical uncorrelated particles, which is recovered in 
the tunneling limit [2] . The experimental determination of 
FCS, unfortunately, is a difficult task. Motivated by the 
expertise developed for current-noise measurements, one 
possibility is to adopt the strategy of building up the FCS 
through the measurement of the various cumulants of the 
distribution, noise being the second one. The third cu- 
mulant was indeed directly measured in Refs. [3,4]. The 
alternative possibility is the direct determination of the 
entire FCS [5]. Recently, various proposals have been put 
forward in this direction making use of the high sensitiv- 
ity to fluctuations of current-biased Josephson junctions 
(JJ) [8-17]. In this Letter we are interested in the sec- 
ond strategy and we analyze the case in which the mea- 
surement is performed applying very short current pulses, 
which allows to access a regime never explored before. 

According to the RCSJ model [18], the dynamics of the 



phase difference across a current-biased JJ is equivalent 
to that of a particle in a washboard potential, the phase 
playing the role of the spatial coordinate. In the harmonic 
approximation, the wells of such a tilted sine potential are 
characterized by the plasma frequency o; p . The wells are 
separated by barriers if the current applied to the JJ is 
smaller than the superconducting critical current I c . In 
such a case, the phase-particle gets trapped in a well of the 
potential, causing the voltage across the JJ to vanish (su- 
percurrcnt state). Escape from the well, which will cause 
the development of a finite voltage across the JJ (transi- 
tion to the resistive state), can occur through two mecha- 
nisms. At low temperatures (T ^ hujp/kB), the only pos- 
sibility is macroscopic quantum tunneling (MQT) through 
the barrier which separates two successive wells. The ef- 
fect of thermal fluctuations on the escape rate of the MQT 
was considered by Martinis and Grabert [19], who found 
an exponential enhancement of the rate proportional to 
T 2 for ohmic damping. The second mechanism, thermal 
activation (TA), is due to excitations of the particle that 
will allow the latter to overcome the barrier top. In the 
absence of perturbations, these are produced by thermal 
fluctuations at large temperatures (T > hajp/ks)- 

In fact, the probability of escape from a potential well 
is sensitive to perturbations affecting the system, such 
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as bias-current fluctuations. This is precisely the phe- 
nomenon current fluctuations detection is based on. In- 
deed, by adding to the bias current, assumed to be con- 
stant, the fluctuating component of the current under in- 
vestigation, it has been shown that the monitoring of the 
voltage appearing across the JJ can be used to characterize 
such fluctuations, either in the TA or MQT regime. Re- 
garding the former, in Ref. [8] an array of overdamped JJ 
which realizes a threshold detector was proposed, whereas 
in Ref. [12] a method was presented for the determination 
of the third cumulant. On the other hand, MQT was ex- 
ploited for the determination of the FCS in Ref. [9] and of 
the fourth cumulant in Ref. [11]. Successful measurements 
in the TA and MQT regimes were very recently reported, 
respectively, in Ref. [13] for the second cumulant, and in 
Rcfs. [15, 16] for the third cumulant. In this paper we are 
interested in hysteretic (underdamped) JJ in the MQT 
regime. Inspired by the experiments [20-22], we consider 
the situation in which the measurement is performed by 
applying a sequence of very short current pulses (of du- 
ration At of the order of a few nanoseconds) [23]. This, 
together with a proper filtering of the frequency compo- 
nents of the current fluctuations, allows to realize a sep- 
aration of time scales, as we shall argue in the following 
(see Fig. 1(a)). 

Current fluctuations whose frequency uj is larger than 
the plasma frequency lu p give rise, even at zero temper- 
ature, to thermal-like activation and represent the non- 
adiabatic regime. For fluctuations of frequency smaller 
than the plasma frequency, escape from the well can only 
occur through MQT, which turns out to be exponen- 
tially sensitive to current fluctuations. This is the adi- 
abatic regime since the particle remains in the ground 
state of the well. Within the adiabatic regime, one can 
realize two different situations depending on whether one 
selects the current fluctuations which occur on a time 
scale longer (low frequency regime) or shorter (high fre- 
quency regime) than At. For large At, of the order of //s, 
only the high frequency (HF) regime is accessible, since 
long-time-scale fluctuations do not contribute appreciably 
(see Refs. [9, 15]). In this case the selection is obtained 
by placing a low-pass filter of bandwidth u>h <C lo v be- 
tween the source of current fluctuations and the J J. In this 
Letter we shall focus on the low-frequency (LF) regime, 
where At is of the order of a few nanoseconds [20-22] and 
< 2w/ At -C w p in order for the HF components not to 
affect the escape probability [24]. The current will fluctu- 
ate for different pulses, remaining constant within a single 
pulse. We shall discuss how the distribution of current 
fluctuations can be extracted from the escape probability 
histogram and, in the small skewness limit, how the for- 
mer depends on the second and third cumulants. We shall 
furthermore show that in some cases the influence of the 
current distribution in the LF regime on the probability 
histogram is more transparent than that relative to the HF 
regime. In addition we discuss in detail two different ex- 
amples. We shall assume that the current fluctuations are 



produced by: i) a tunnel junction; ii) a source of random 
telegraph noise [25]. 

The system. — The system we consider (see 
Fig. 1(b)) is a dc-SQUID, composed of a superconducting 
loop with two identical J J. We consider a loop of small 
inductance, such that the SQUID is equivalent to a sin- 
gle JJ with a flux-tunable critical current. The SQUID 
is biased by a constant current source I and coupled to a 
voltage amplifier. The fluctuating component 6I(t) of the 
current originating from a mesoscopic conductor is fed into 
the SQUID through a proper filtering circuit (not shown 
in the figure) [26]. Such a circuit represents the environ- 
ment of the SQUID which must be taken into account in 
the interpretation of the results of an experiment. Indeed, 
the environment will introduce additional contributions to 
the current fluctuations under investigation even exerting 
a back-action on the system [28,29]. Since the exact effect 
is determined by the details of the circuit, in the following 
we shall disregard the effects of the environment and focus 
only on the mere effects of current fluctuations. 

Current fluctuations arc characterized by the distri- 
bution p(SI) which can be written as a Fourier trans- 
form p(SI) = dke~ lkSI <f>si{k) of the characteristic 

moment-generating function <j>si(k) = cxp(^^L 2 c n ), 
where c„ are the cumulants of the distribution. In partic- 
ular, c-2 is related to the width (c2 = (SI 2 )) and C3 to the 
asymmetry of the distribution (03 = (5 1 3 )). 

The J J, which is assumed to be underdamped, is oper- 
ated as follows. A sequence of current pulses, of duration 
At and amplitude /, is applied to the JJ, and the voltage 
across the JJ is monitored. If, during the time interval 
At, a finite voltage across the JJ develops, the event is 
recorded. The ratio between the number of transitions and 
the total number of pulses applied is defined as the escape 
probability, which, as a function of the average amplitude 
/ of the bias current, gives rise to a histogram P(I) (see 
Fig. 1(c)). Such an escape probability histogram is expo- 
nentially sensitive to perturbations affecting the system, 
such as current fluctuations SI(t) produced by the meso- 
scopic element or fluctuations of the magnetic flux thread- 
ing the SQUID. We assume that temperature is small such 
that we can neglect thermal fluctuations. 

For very short pulses we access the LF regime, in which 
the current is constant within a single pulse, but fluctu- 
ates for different pulses. Assuming ergodicity, the escape 
probability is given by averaging the probability in the 
absence of fluctuations Pq over the distribution of current 
fluctuations p(SI): 

/+00 
d(6I)p(SI)P Q (I + 61), (1) 
-00 

where Pq(I) = 1 - e~ rA *, T = T(I,<&) being the escape 
rate in the MQT regime given by T(I) = A(I) exp(-B(I)), 

where A(I) = U^g^ and B(I) = 

in the limit of low dissipation [30]. The plasma ffe- 
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Fig. 1: a) The frequency arrow: fluctuations with frequency 
to <C Up produce an adiabatic perturbation of the washboard 
potential, thus affecting the escape probability with exponen- 
tial sensitivity. High frequency (HF) fluctuations (u) > 2tt/ At), 
occur over the duration of a pulse giving rise to an averaged es- 
cape rate. On the contrary, for low frequency (LF) fluctuations 
(uj < 2n/ At), the current does not fluctuate during a pulse, but 
changes from pulse to pulse giving rise to an averaged escape 
probability. In the LF regime we assume < 2n/At <C u) p , 
where ub is the bandwidth. Non-adiabatic perturbations are 
due to fluctuations with frequency uj > u p which give rise to 
thermal-like activation, b) The system consists of a current- 
biased dc-SQUID connected to a voltage amplifier. 51 repre- 
sents the fluctuating component of the current flowing through 
a mesoscopic element, whose current fluctuations distribution 
is under investigation, c) Escape probability histogram ob- 
tained by varying the bias current. 
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where 
/2Cj is the 



quency is given by haj p = \/8EJEq 

Ej = hl c /2e is the Josephson energy, Ec = e 
charging energy and 7 C is the critical current of the SQUID 
for zero magnetic flux. AU denotes the barrier height and 

it is defined as AU = 2 



E, 



3/2 



The distribution p(5I) can be derived directly by de- 
convolving the measured escape probability, once Pq is 
known. By Fourier transforming Eq. (1) one obtains the 
characteristic function <fi(k) = P(k)/Po(k) where P (Pq) 
is the Fourier transform of P (Pq). The various cumu- 
lants of the distribution can then be calculated through 
the derivatives c„ = (— i) n Jp- log</>(fc)| fe _ . In the case 
where Pq has a very sharp transition [Po(I) — 0(1 — Iq)] 
the distribution can be simply obtained by a derivation 



p(x) 



dP(I - x) 
dx 



(2) 



so that the central moments of the distribution arc ob- 
tained = J dx x n p(x). We wish to remark that 
in the LF regime the low frequency component only will 
contribute to the moments and hence to p. For example: 
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x(A/(t 1 )A7(t 2 )A7(0)). 



Higher frequency (adiabatic) components are filtered out 
by the external circuit. 



Small skewness limit. 

higher moments of current, 
skewness 7 of the current 

c 3/ c 2^ 2 *^ I); one obtains 



Neglecting the fourth and 
and in the limit where the 
distribution is small (7 = 



p(5I) 



1 



V27TC2 



2c 2 



(7) 



We shall refer to the escape current 7 esc as the current rel- 
ative to P(7 esc ) = 1/2, and to the width of the transition 
as A7 = 7o.9 — 7o.i, with I x defined by P(I X ) = x. 

For the sake of definiteness, we shall assume that Pq 
has a very sharp transition (with respect to the values 
of c 2 /7o considered), so that we approximate the proba- 
bility with a theta function Po(I) = 8(1 — Iq). Let us 
now discuss the behavior of the probability histogram in 
the LF regime. For C3 = (Gaussian noise), the escape 
current does not depend on c 2 , while the width of the 
transition A7 increases linearly with y/c2. More precisely, 
A7 = 2V2liErf _1 (4/5), where Erf" 1 indicates the inverse 
of the error function. Now, a finite skewness is expected 
to produce an increase (decrease) of the escape current for 
7 > (7 < 0). We calculate the escape probability and 
determine 7 esc and A7 as a function of 7 for different val- 
ues of c 2 . As shown in Fig. 2, 7 esc increases linearly with 7, 
the slope being proportional to •J^./Iq in such a way that 
(7 esc — 7o)/y / c 2 = const • 7. Since A7 at 7 = depends on 
the value of c 2 , in Fig. 3 we report the plot of the relative 
width change A7 r = A7 — A7(7 = 0) as a function of 7. A 
quadratic dependence on 7 is found, with a coefficient pro- 
portional to y/c^/Io so that A7 r / v / c 2 = const • 7 2 . Note 
however that the dependence is very weak: 0.01% escape 
current change and 0.1% transition width change in the 
range of 7 considered. We wish to mention that the results 
in the case of a finite-width transition, in the absence of 
fluctuations, arc qualitatively equal. For example, in the 
experimentally relevant case in which A7i nt /7o = 5 x 10 -3 , 
where A7j nt is the intrinsic width of Pq, wc find the same 
dependence of 7 esc and A7 on 7. The only quantitative 
difference is that the values of A7 are shifted upwards, 
while the curve of 7 osc is slightly shifted to lower values. 

It is interesting to compare the above results with the 
ones relative to the HF regime, in which the current fluc- 
tuates within a single pulse. In such a case the escape 
probability is given by averaging the escape rate over the 
noise distribution p(SI) of Eq. (7): 



P(7) = 1 



= 1 P -T(/))At _ 



1- 



,~ /- °° d(SI)p(SI)r(I+SI)At 



(8) 



F. Taddci and F. W. J. Hekking 




Fig. 2: I eS c/Io as a function of 7 for several values of C2. A Fig. 3: The normalized relative width change A7 r /7o is plotted 
linear dependence is found with slope proportional to ^02/ Iq. as a function of 7 for several values of C2. 



It turns out that the width of the transition AI does not 
depend on 7 so that the skewness affects the escape current 
only. However, since the influence of C2 on the escape 
current is strong, it is difficult to isolate the shift in the 
escape current due to a change of 7. 

In Ref. [27] it was shown that, in the HF regime, the ef- 
fect of the third moment on the escape rate can be singled 
out by inverting the polarity of the current /. Such in- 
version results in a mere shift of the escape current which 
allows, in principle, to evaluate the value of 7. We apply 
a similar analysis to the LF regime finding, interestingly, 
that it is the escape probability as a whole that is affected 
by the presence of the third moment. For the average 
escape probability we obtain 

1 f°° 

P ± (I) = -= / dxP {I + ^x) exp(-x 2 /2) T 



7 I 



2 J2n 



dxP (I + y/c^x)x{\ - x 2 /3) exp(-a; 2 /2), (9) 



where the subscript ± refers to positive (negative) current 
polarities. The values of 7 and C2 can be determined if Pq 
is known. In particular, when Pq can be approximated by 
a theta function, the value of 7 can be determined from 



p+(i)-p-(i) = -r 



2c 2 V2ir 
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(10) 



while C2 can be determined from the sum P + + P~ . It 
turns out that the effect of the inversion of the current po- 
larity is smaller in the LF regime than in the HF regime. 
However, while for LF 7 can be evaluated directly by com- 
paring two histograms (see Eq. (10)), for HF 7 can be 
determined only indirectly from the escape rates and the 
histogram derivative [27]. Finally we note that spurious 
current asymmetries may arise in the SQUID as a result of 
parasitic magnetic fields. To avoid this phenomenon one 
should replace the SQUID with a JJ as detector. 

Tunnel junction. — Let us consider a concrete exam- 
ple in which the mesoscopic element is a tunnel junction. 



The current flowing through it is distributed according to 



/Otun(-f) 



e r (ii+i) 



(11) 



where N = T/el is the mean number of electron trans- 
ferred in a time T, I is the mean current, e the electronic 
charge, and T is the Euler Gamma function. Despite the 
fact that the FCS of a tunnel junction is Poissonian, with 
all cumulants equal to N, the cumulants of the current 
distribution p tun are given by c n = J(e/T)™ _1 , i.e. de- 
creasing with increasing n. It is interesting that for N as 
large as 100, ptun is easily distinguishable from the Gaus- 
sian distribution of equal C2 . Such small values of N can 
be achieved, in our detection scheme, since the time T can 
be identified with 2ir/uj-t,, which is taken to be of the order 
of the pulse duration At. What limits the smallness of TV 
is, firstly, the minimum achievable value At, and secon- 
darily the maximum value of lu^, which must be smaller 
than uip. As an example, in Fig. 4, we plot the function 
Ptun(5I) = ptun(SI + I) {i. e. the distribution shifted by 
mean current I) together with the corresponding Gaussian 
(red curve), for At = 1 ns and I = 10 nA, so that N = 65 
and 7 = N^ 1 / 2 ~ 0.12. For clarity, in the inset we plot 
the same curves on logarithmic scale. 

The simplest way to measure ptun in the LF regime is 
through Eq. (2), valid in the limit of a step- like transition 
in the escape histogram Pq. The latter condition requires 
that the width of the transition in the absence of noise 
Ai nt is much smaller than the width of the distribution 
y/c~2. However, this is hardly possible with realistic pa- 
rameters, since the above two quantities arc typically of 
the same order. Indeed, in a typical experiment Ai nt /7 C is 
of the order of 1 x 10 -3 , whereas with / = 10 nA, I c = 1 
HA and At = 5 ns one finds ^/I c = 6 x 10~ 4 [31]. Sur- 
prisingly, the derivative of the probability histogram Ptun, 
obtained by Eq. (1) using ptum turns out to be distinguish- 
able from the corresponding derivative of ^baussj relative 
to a Gaussian distribution of same width. This can be 
appreciated in the left-hand-side of the plot reported in 
Fig. 5, where the above two quantities are plotted as a 
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Fig. 4: Distribution of current fluctuations ptun{51) relative to 
a tunnel junction with At = 1 ns and I = 10 nA (black curve) 
and of the corresponding Gaussian of equal width (red curve). 
In the inset we plot the same curves in logarithmic scale. 



function of the normalized current fluctuations. Even in 
this far-from-ideal situation, the effect of higher cumulants 
is weak but present [32]. It is worthwhile mentioning that, 
in order to make use of Eq. (2), one needs to amplify the 
cumulants by a factor a, for example through an external 
circuit. For the figures reported above, a value of a = 100 
appears to be sufficient for a determination of ptun- We 
shall sec in the next paragraph that the LF regime is much 
more suited for peculiar distributions. 

Flux random telegraph noise. — Next we consider 
the situation where the flux threading the dc-SQUID in 
Fig. 1 is affected by random telegraph noise (RTN) and, for 
defmiteness, where the current is noiseless. The possibility 
of such a source of noise has been considered in Ref. [33] . 
The time-dependent flux can be written as </>(t) = <j> + 
A(f> cos [tt ^2 k 9(t — tk)] , where <j> is the average flux, A<^ is 
the displacement and are randomly distributed times. 
The distribution of flux fluctuations 5<fi is given by 

Paux(<ty) = \ W4> - A0) + 5(5cf> + A<j>)\ , (12) 

with non-zero even central moments {({8<j)) 2n ) = (A<fi) 2n ) 
and characteristic function x(A) = i (e iAA< ^ + e~ iAA< ^) . 
The time-dependent correlation function is given by 
(5cj)(t)6(f){0)) = (A0) 2 e - 2 l*IP, where p is the probability of 
transition between the two states (+A0/2, — A0/2) per 
unit time. The noise power (double Fourier transform) is 

therefore a Lorenzian S(uj) = 4p2+u 2 ■ ^ n ^ s case 

one can choose the values of p and At such that flux fluc- 
tuations has only a LF component. Indeed, if we assume 
that 1/At 3> Pi the escape probability is given by 



d(<tyVflux(<ty)Po(/,0 + <ty), (13) 



in the case where there is no current noise. In Eq. (13), the 
flux dependence in Po(I,(j)) is introduced by taking into 

account the fact that I c = 2/ c j,j cos (j^fj , where 7 c jj 

is the critical current of a single JJ. In Fig. 6, P(I, 4>) is 
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Fig. 5: Logarithmic plot of the derivative of P tun relative to a 
tunnel junction (black curve) for realistic parameters describ- 
ing the J J (see text). In red is plotted the corresponding curve 
for a Gaussian current distribution of equal width, which can 
be distinguished from the black curve, in the range of values 
considered, at least for negative fluctuations. 



plotted for <f>/<j)a = 1/5, A(j}/c/) = 1/60, I c = 1/iA, At = 5 
ns and Eq/Ej = 10~ 4 (<fo = h/2e being the flux quan- 
tum). The red curve refers to the absence of noise, the 
black curve to RTN and blue curve to Gaussian noise (of 
equal width). The escape histogram relative to the RTN 
presents a double step structure, well distinguishable from 
the histogram relative to Gaussian noise. The double step 
histogram can be understood as due to averaging two his- 
tograms with different escape current, once a finite average 
flux cj> is assumed. However, the effect of RTN noise in the 
HF regime leads to a decreased escape current. The gen- 
eral feature of the HF regime is confirmed also in this case: 
the effect of higher cumulants on the histogram reduces to 
a shift in the escape probability, see also Ref. [11]. 

Conclusions. — To summarize, we have shown that 
the escape probability histogram of an underdamped 
current-biased dc-SQUID, operated with very short cur- 
rent pulses, can be used for retrieving information on the 
distribution of current fluctuations which characterizes a 
given conductor. The latter is connected to the dc-SQUID 
through some filtering circuit, which we do not consider, 
for the sake of defmiteness, but would affect the measure- 
ment. More precisely, this method is sensitive to the low 
frequency component of the fluctuations, namely up to 
1/At (At being the duration of the pulse, of the order 
of a few nanoseconds). We have first analyzed the case 
in which the distribution presents a small third moment 
(small skewness limit), while all other moments are van- 
ishing. In particular, we have addressed the functional 
dependence of escape current and width of the transition 
with respect to the skewness, and found that, even for 
short pulses, the value of the skewness can be in principle 
evaluated by inverting the polarity of the bias current. We 
have then analyzed two concrete examples: in the first one 
we have assumed that the current fluctuations originates 
from a tunnel junction and we have evaluated the pos- 
sibility of determining the whole current distribution, by 
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Fig. 6: Probability histograms with flux LF noise with 4>/4>o — 
1/5, A0/<j!.o = 1/60, 7 C = lfiA, At = 5 ns and E c /E 3 = 10~ 4 . 
Black curve is relative to RTN, blue curve to Gaussian noise 
and red curve to no noise. 7 is expressed in units of 27 c jj. 

simulating the probability histogram under realistic con- 
ditions. We wish to mention that the same analysis can be 
performed for any mesoscopic element. In the second ex- 
ample, we have considered flux fluctuations characterized 
by a random telegraph noise distribution, finding that pe- 
culiar distributions such as the latter produce a very clear 
signature, not present if the measurement is performed 
with long current pulses. Finally we wish to stress that 
we focused the analysis to current pulses only, even though 
it remains valid for the case of flux pulses. 

* * * 
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